BONN-TH-2010-01 



Direct Integration and Non-Perturbative Effects 

in Matrix Models 

o 

', Albrecht Klemm a |l|, Marcos Marino 6 ]^ and Marco Rauch a jf| 

P-h 

a Bethe Center for Theoretical Physics, 
Universitdt Bonn, Nufiallee 12, D- 531 15 Bonn, Germany 



b D epartement de Physique Theorique et Section de Mathematiques, 
Universite de Geneve, CH-1211 Geneve, Switzerland 



Abstract 

qq ! We show how direct integration can be used to solve the closed amplitudes of multi-cut 

matrix models with polynomial potentials. In the case of the cubic matrix model, we give 

CN ■ 

O explicit expressions for the ring of non-holomorphic modular objects that are needed to 

express all closed matrix model amplitudes. This allows us to integrate the holomorphic 
anomaly equation up to holomorphic modular terms that we fix by the gap condition up 

•rH . 

^ . to genus four. There is an one-dimensional submanifold of the moduli space in which 
the spectral curve becomes the Seiberg-Witten curve and the ring reduces to the non- 
holomorphic modular ring of the group T(2). On that submanifold, the gap conditions 
completely fix the holomorphic ambiguity and the model can be solved explicitly to very 
high genus. We use these results to make precision tests of the connection between the 
large order behavior of the 1/N expansion and non-perturbative effects due to instantons. 
Finally, we argue that a full understanding of the large genus asymptotics in the multi-cut 
case requires a new class of non-perturbative sectors in the matrix model. 
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1 Introduction and Results 



In this paper we propose direct integration as a new method to solve the closed amplitudes 
for multi-cut matrix models with polynomial potentials. More precisely, we calculate the 
closed partition function of such matrix models 



perturbatively in the genus g, but exactly in the 't Hooft parameters S_. Exact means that 
the F g (S_) are given in terms of period integrals of the spectral curve and can be written 
explicitly in terms of modular forms of subgroups of SP(2g, Z). 

Direct integration refers to a method of solving the holomorphic anomaly equation [B] 
using the modular transformation properties of the amplitudes under the monodromy 
group of the spectral curve. This method has been developed in the context of topo- 
logical string theory in [5U E21 I2S1 El [23 I2E]- The fact that the holomorphic anomaly 
equations govern such matrix models was suggested by the large N duality of [TJ5]. In 
this duality, type B topological string amplitudes on certain local Calabi-Yau spaces turn 
out to be encoded in the 1/N expansion of matrix model partition functions. Therefore, 
the holomorphic anomaly of the topological string naturally carries over to these matrix 
models as first pointed out in [30]. It has been shown much more generally in [22] that 
the holomorphic anomaly equation is valid for all matrix models which are solvable by 
the method of [23]. 

The holomorphic anomaly equation relates anti-holomorphic derivatives of the closed 
amplitudes F g (S_) at genus g to lower genus amplitudes F h<g (S_), in a recursive way. Since 
only the anti-holomorphic derivative is specified by the equations, the procedure leaves a 
holomorphic ambiguity, i.e. F g (S_) = F g h (S_) + f g (S_) splits into a non-holomorphic term 
F g h (S_), which is determined by the holomorphic anomaly equation, and the holomorphic 
ambiguity f g (S_), which must be fixed genus by genus by using modular properties and 
boundary conditions at special points in the moduli space. The modular transformation 
properties imply that the amplitudes are generated by a finite ring of modular forms, 
which have holomorphic as well as non-holomorphic generators. Modularity and the 
holomorphic anomaly equation imply that the total amplitude F g (S_) is a polynomial in 
these generators whose degree grows linearly with the genus. The ambiguity f g (S_) is a 
polynomial generated by the smaller ring of holomorphic generators. The finite number 
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of coefficients in this polynomial must be fixed by boundary conditions. 

In this paper we find that the gap conditions, which where investigated in non- 
compact [30j EHJ HI El] and compact Calabi-Yau backgrounds [32} [26J [271 [28], provide 
enough independent boundary conditions to fix the ambiguity (and hence the amplitudes) 
completely. Following [2H] we refer to this property as integrability of the holomorphic 
anomaly equation. 

The large N duality relating matrix models and topological strings gives a natural 
geometric interpretation to the algebraic objects describing the planar limit of the matrix 
model [19]. The spectral curve y(x) of the matrix model (which, in the case of polynomial 
potentials, is a hyperelliptic curve) describes the distribution of eigenvalues in the planar 
limit, and in the topological string dual it describes the nontrivial part of the Calabi-Yau 
geometry. We derive the modular ring starting from the Picard-Fuchs equations governing 
the periods of the form Q = y(x)dx. This is a general method^), and since we expect that 
the gap boundary conditions fix the ambiguity, our approach should apply to general 
multi-cut matrix models with polynomial potential. 

Of course, the formalism of [23] gives in principle all the genus g free energies of generic 
multi-cut matrix model in terms of universal formulae on the spectral curve. The price 
to pay for such a general approach is that its detailed implementation is in practice very 
involved. Even in two-cut models, going beyond genus two with the methods of [23] is 
not very feasible. In contrast, direct integration becomes very powerful when the spectral 
curve and its modular group are simple. 

In this paper, in order to illustrate the method of direct integration, we focus on the 
two-cut matrix model with a cubic potential. In this model the N eigenvalues split in two 
sets N = N\ + N2 and condense in sets near the two critical points of the potential. This 
leads to the cuts in the spectral curve shown in figure [TJ There are two independent 't 
Hooft couplings Si = g s Ni, i = 1,2, which correspond to the integrals of Q over the two 
cuts. As shown in [19], the planar free energy of this matrix model, F (S±, S2), calculates 
the exact superpotential W e g of an J\f — 2 U (M) supersymmetric gauge theory broken 
down to an M = 1 gauge theory U(M\) x U(M 2 ), by a cubic three-level superpotential 
in the adjoint [19] (notice that iV» are unrelated to Mi). The higher genus amplitudes 
F g (Si, S2) in the matrix model arise as generalized couplings in a non- commutative de- 

4 For example, a meromorphic modular form of weight k of SL(2,Z) or a congruence subgroup fulfills 
a linear differential equation of order k + 1 in the total modular invariant [53] . 
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formation of the M = 1 gauge theory [44] . 

Certain aspects of the original M = 2 theory can be recovered from the J\f = 1 theory 
by breaking the gauge symmetry to the Cartan subgroup and taking the limit in which 
the superpotential vanishes [13]. When the gauge group is 577(2), a cubic superpotential 
is enough to go to the Coulomb branch. This implies that various quantities appearing in 
the Seiberg-Witten solution of pure M = 2 super Yang-Mills theory [47] can be obtained 
from a matrix model calculation with a cubic potential, and on the slice Si = — S2. These 
include the gauge coupling [T7] and the R\ gravitational coupling [3U [18]. In fact, the 
spectral curve of the cubic matrix model on that slice is identical to the Seiberg-Witten 
curve [17] . Since the modular group of this curve is particularly simple, direct integration 
becomes an extremely powerful method to calculate the F g (Si, —Si), as we show in section 
® 

On the other hand, in the SU(2), M = 2 gauge theory there is an infinite number of 
couplings F g (a), g > 2, which describe the gauge-gravity couplings F^~ 2 R 2 + involving the 
graviphoton field strength F + . These couplings appear naturally in Nekrasov's partition 
function [42] and they can be also obtained by using the holomorphic anomaly equations. 
This was shown for the pure gauge theory and SU(2) with matter in [301 126] and [3Tj 
respectively. However, it was noticed in [M] that these higher genus couplings F g (a) do 
not agree with the higher genus F g (Si, S2) obtained in the cubic matrix model and then 
restricted to the slice S2 = —Si. This disagreement is due to the fact that the Seiberg- 
Witten differential Asw differs from the natural differential Q on the spectral curve of 
the matrix model. In contrast, r and F\ only depend on the spectral curve, and not on 
the differential, and therefore are the same in both cases. In j35j 09] matrix models are 
derived which encode all M = 2 gauge theory amplitudes F g for arbitrary g, however one 
has to introduce potentials involving polylogarithms and their quantum generalizations. 

An interesting application of our computation of the couplings F g (Si,S2) at high 
genus is the study of non-perturbative effects in matrix models and their connection to 
the large order behavior of the 1/7V expansion. It is well-known that, in many quantum 
systems, there is a connection between perturbation theory at large orders and instantons 
(see for example [3?]). In matrix models, instanton configurations correspond to the 
tunneling of eigenvalues between different saddle points p~41 148] . A detailed analysis of 
these configurations for off-critical, one-cut matrix models can be found in [40], which 
verified the connection to the large order behavior of the 1/7V expansion in detail in 
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some nontrivial examples. In this paper we explore this connection in the two-cut matrix 
model. On the one hand, we find that the large order behavior is controlled at leading 
order by the action of a single eigenvalue tunneling from one saddle-point to the other, in 
agreement with the general ideas put forward in [TU HH HQ]. On the the other hand, we 
argue that a full understanding of this connection requires new non-perturbative sectors 
which have not been yet identified in the matrix model. The existence of these sectors is 
also suggested by a recent analysis of the asymptotic behavior of the instanton solutions 
of the Painleve I equation [25]. We conjecture that these sectors might involve topological 
brane-antibrane systems. 

The paper is organized as follows. In section [2] we give the general ideas of the direct 
integration method and of the modular covariant approach of [I]. In section 3 we describe 
in detail the geometry underlying the two-cut cubic matrix model, we set up the direct 
integration formalism and we analyze the boundary conditions. In section 4 we point 
out that the modular covariant formulation is most powerful on the one- dimensional slice 
of moduli space S\ = —5*2, and we develop direct integration on this submanifold in 
terms of non-holomorphic modular forms of I\2). Section 5 is devoted to the analysis of 
instanton effects in the two-cubic matrix model, and their connection to the large order 
of the 1/N expansion. Finally, the Appendices collect some useful information about the 
cubic matrix model as well as results on modular forms which are used in the paper. 



2 Holomorphic anomaly and direct integration 

Below we review very briefly the generic aspects of the techniques of direct integration of 
the holomorphic anomaly equation of [6] 

B- t F g = i<3f (^DjDkFg.! + DjF g - r D k F^j , (g > 1) , (2.1) 

which was derived for Calabi-Yau three- folds in [5H E21 ESI El E3 ES] • In our application 
to the spectral curve £ of a matrix model, the Di are covariant derivatives Di with respect 
to the metric G on the moduli space of the Riemann surface S. 

We note that Cf = C, fk G Jj G kk , where C ijk can be derived from the holomorphic pre- 
potential F as C^k = DiDjdkF Q . The prepotential F , the metric G^ and flat coordinates 
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S_ can all be derived from the period integrals 

ft, / ft), i = !,...,<?(£) (2.2) 



/ 

over a symplectic basis (a\ 6j) of i/^X^Z). In particular, given a point in the moduli 
space, one can make a choice of this symplectic basis, so that suitable flat coordinates are 
defined by 

S l = /ft (2.3) 



a' 



while the hi periods IT fulfill 

n=^ (2 4) 

These relations determine the prepotential F up to an irrelevant constant. We define the 
r matrix of the Riemann surface as 

r--^L- (2 5) 

u "as<asr 1 ] 

The matrix Im(r)jj is positive definite, and it gives the metric on the moduli space of the 
model. Equivalently, the metric can be obtained form the Kahler potential 

K = ^-(n i S T -n,S>). (2.6) 

On Riemann surfaces the period integrals can often be directly performed. Alterna- 
tively it might be useful to derive the Picard-Fuchs equations and reconstruct the periods 
as linear combinations of their solutions. Much of the above has been spelled out in the 
context of the Riemann surfaces for the B-model of topological string theory on non- 
compact Calabi-Yau in [29]. The relevant compact part of the geometry is given by a 
Riemann surface and a meromorphic differential, which comes from reducing the holo- 
morphic (3, 0)-form on the Riemann surface. After identification of the former with the 
spectral curve £ and the later with the form ft, we can use the formalism discussed in [29] . 

One property of the matrix model geometry is that the periods over the a-cycles do 
not fulfill the relation YH=i & l = 0- Usually this relation is inherited by the periods of 
holomorphic forms due to the homological relation of the cycles. However, in matrix 
models one has Yll=i & l ^ N, because ft has one non-vanishing residue outside the cuts. 
This leads to one algebraic relation between the periods in terms of the r parameters, 
which for the r = 2 case (the cubic matrix model) is expressed in eq. (13 .28 p . The property 
of a non-vanishing residue is shared with Seiberg-Witten theories with matter [31] and 
certain non-compact Calabi-Yau geometries with more than one Kahler class j29j H] . 
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2.1 Direct integration 

The so-called propagator plays a decisive role in the solutions of the B-model [B] . For the 
formalism on the Riemann surface £ one needs only one typd^l of propagator defined 
by 

where i,j = 1, . . . , r and r is the number of parameters in the model. Following [6] it can 
be shown that the F g can be written as 

39-3 

F 9=J2 (S) S^ 2 . . . S^M'm (2.8) 
|/|=o 

where the f g j{S_) are holomorphic tensors of the moduli. The most important property 
of the S ij is that 

dF 

d i F g = Cf^. (2.9) 

If one assumes linear independence of the as functions of S_, it follows from this 
property that (12. ip can be rewritten as a set of equations 



dF^ 

dSv 



\ [DjDkFg^ + J2 D j F 9-rDkF r ) , (<? > 1) . (2.10) 



r=l 



These equations can be integrated algebraically, provided that the r.h.s. can be expressed 
in terms of the S 1 ^ contracted by holomorphic tensors as in the r.h.s of ( 12.81) . This is 
possible since the following closing relations are fulfilled due to special geometry [6j [3] 



D^ 1 = -C mm S km S ln + ft, (2.11) 
r| = -C^S" 1 + ft , (2.12) 



1 



d i F 1 = -C ijk S Jk + Ai. (2.13) 

Here the /• ■ , ff l and Ai are holomorphic ambiguities, which must have the same trans- 
formation properties as the expressions on the left-hand side. These ambiguities are due 
to the fact that (j2.7p defines only up to an holomorphic tensor. Different choices are 
possible and lead to a redefinition of the f g j in (I2.8p . As we mentioned above the periods 



3 In the threefold cases there are three types , S l and S. 
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are not algebraically independent, see for example (13.28)) . As a consequence it is possible 
to make a choice for the above ambiguities so that for a given i one has S lk = 0, Wk, 
i.e. the matrix of propagators has effectively only rank p = r — 1. We call the auxiliary 
parameter t . There may be more auxiliary parameters stemming from the independent 
non- vanishing residua of Q. If there are k such residua, the rank is reduced to p = r — k. 

Whether one works with the redundant or the reduced set of propagators the equation 
(12.101) can easily be integrated w.r.t. S 1 *- 7 and F g becomes of degree 3g — 3 in the . This 
is an efficient way to solve the recursion, but at each step one still has to determine the 
holomorphic ambiguity. 

2.2 Modular covariant formulation 

It is possible to relate the non-redundant set of propagators to quasimodular forms. In 
particular, in the holomorphic polarisation, the following properties derived in [1] hold: 

1. F g (S_) is invariant under the monodromy group T of the Riemann surface E. 

2. F g (S_) is an almost- holomorphic modular function, i.e. its non-holomorphic depen- 
dence is encoded solely in ((r — f )~ 1 ) /J , where I, J — 1, . . . , p and tjj is the standard 
matrix valued modular parameter living in the Siegel upper half spaced 

3. The non-holomorphic dependence combines always with quasimodular forms E IJ to 
give almost-holomorphic modular forms 

E IJ = E IJ (T) + ((r-f)- 1 ) IJ . (2.14) 

Here we defined E IJ (r) as derivative of ^ji^i(r). The anomaly equation of [6] for 
.Fi(t) implies that this is a non-holomorphic modular invariant 



Fi = - log 



det*(Im(Tij)) ($ fc (f)$ fc (r)) a j (2.15) 



under the monodromy group V. $fc(r) is a holomorphic Siegel modular cusp form of 
weight k which vanishes at the discriminant A of the Riemann surface. It transforms 
as $fc(r 7 ) = det (Cr + D) k § k (r), where r 7 is given by 

t j = (At + B)(Ct + D)~ 1 , 7 =M *) e S P (2,,Z). (2.16) 



= r — k, where k is the number of independent non- vanishing residua of fi. 
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Such modular forms exist for all genus and can be written as products of even theta 
functions [36]. The exponent a will make the argument of the log invariant and 
the vanishing order at the discriminant ^log(A). For an elliptic curve $fc(r) is 
typically the Dedekind ^-function. However, if the subgroup V allows for several 
cusp forms, $fc(r) can be a suitable multiplicative combination of them. In virtue 
of the definition E IJ transforms as a Siegel modular form 

E ij (t 7 ) = (Ct + DY K (Cr + D){E K \r) . (2.17) 
4. F g (S_) can be expanded as 

39-3 

F 9 = ^ Uh,..ME Ilh ■ ■ ■ E 1 "-^ . (2.18) 

|J|=0 

Note that f g j has to compensate for the modular transformation of r and can in principle 
be expressed through holomorphic modular forms. 



3 The two-cut cubic matrix model 

As shown by Dijkgraaf and Vafa in [19], the B-model topological string theory on certain 
non-compact Calabi-Yau geometries is captured by a matrix model. The matrix model 
is the n-cut matrix model with potential W(x), while the Calabi-Yau geometry is the 
following hyper surf ace in C 4 

uv = y 2 - [W\xf + f{x)). (3.1) 

Here, f(x) is a polynomial of degree n— 1 that splits the n double zeroes of W'(x) 2 , see [38] 
for a detailed review. In the following we will combine the Dijkgraaf- Vafa correspondence 
with known results about the holomorphic anomaly equation in order to give a recursive 
solution of multi-cut matrix models. 



3.1 The geometrical setup 

In this paper we will be interested in multi-cut, Hermitian matrix models. The partition 
function is defined by 

Z = \ - f dMe-i w ( M > (3.2) 
vol(C/(A0) J V ' 
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where W(M) is a polynomial of degree d = n+ 1 in the NxN matrix M. The most general 
saddle point of this model at large N is a multi-cut solution, in which the eigenvalues of 
M condense along cuts 

[a~,af] C C, i — l,---,n, (3.3) 

in the complex plane. The cuts are centered around the n critical points of W(x). One 
way of encoding the planar solution of the matrix model is through its resolvent 

The planar limit of this correlator, denoted by ojq(x), has the structure (see for example 

«*>(?) = ^{W(x)- V (z)), (3.5) 

where 

* = 9sN (3.6) 

is the total 't Hooft parameter, and 

n 

y\x) = (W'(x)f + f(x) = c JJ(x - a~)(x - a+) (3.7) 

i=l 

is called the spectral curve of the multi-cut matrix model. In the matrix model literature 
it is customary to write it as 

y\x) = M{x)^{x), (3.8) 

where cr(x) is a polynomial in x, 

2s 

<r{x)=Y[(x-Xi), (3.9) 

i=l 

and s <n. Of course, if all the roots in (13. 7p are different, s = n and M(x) is a constant. 

Through the large duality of [TJJ], the spectral curve (13.71) describes the nontriv- 
ial part of the target geometry (13.11) . The positions of the endpoints are fixed by the 
asymptotic condition 

u (x) ~ — , x — > oo, (3.10) 

x 

and by the requirement that there are iVj eigenvalues in each cut, 
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In this equation, d is a contour encircling the cut [a { , af] counterclockwise. The partial 
't Hooft parameters Si are defined by 

Si=g s Ni. (3.12) 

Notice that 

n 

In the following we consider a cubic matrix model with potential W given by 

W{x) = ^x 2 + fx 3 . (3.14) 

Since this model has two critical points x = ax, x = a 2 , the generic saddle will be a two- 
cut matrix model. If we write the matrix integral (13. 2p in terms of eigenvalues, we have 
to distinguish two different sets {/ii}i=i,... ^ , {^j}i=i,-,iVa) which are expanded around ax, 
a,2, respectively, and we obtain 

■ 2 ' i=l j=l i<j i,j 

(3.15) 

Since 

W'(a;) = mx + gx 2 = gx^x H J = <7(x — ai)(x — 02), (3.16) 

iy'(x) 2 has two double zeroes at x = ax, a 2 , that are split by the degree one polynomial 

/(z) = Ax + /x (3.17) 

into four roots af, af. Hence, the curve for the geometry/matrix model is given by 

y 2 = W\xf + f = g 2 (x - al){x - af)(x - a 2 )(x - af). (3.18) 

We choose the branch cuts to be along the intervals (a^,af) and (a 2 ,af), cf. Fig. [TJ 
It follows from (13. lip that the 't Hooft parameters for this curve are the periods of the 
one-form 

tt = y(x)dx (3.19) 
around the branch cuts. Following the notation of [12], we have 



1 r a t 1 

Si = — I n, Ui = — I a (3.20) 



27ri ./„- 2ii\ 
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These 't Hooft parameters are functions of the couplings in the potential m, g, and of 
the variables A, fi. Equivalently, they are functions of the branch points af of the quartic 
curve (13.181) . It is convenient to define new variables given by 

z\ = ^(>2 -^i) 2 , z 2 = ^(x 4 - x 3 ) 2 , 

Tfi 

Q = ^{xi + x 2 + x 3 + x 4 ) = -—, (3.21) 
I 2 = \ [(x 3 + X A ) - (X! + x 2 )} 2 = (^j - 2( Zl + z 2 ), 

where we label the cuts more conveniently as 

, af , a^~, af) = (xi, x 2 , x 3 , X4) (3.22) 

and we also have 

4 

a(x)=Y[(x- Xi ). (3.23) 

i=l 

We will use this in order to expand all four periods in powers of Z\ and z 2 . Notice 
that Zi are coordinates that parameterize the complex structure deformations of the local 
Calabi-Yau geometry ( 13. ip . 

Let us consider S\. For this we change variables to y = x — \{x\ + x 2 ) and the integral 
becomes 

g rvA 

Si = — V(y- 2/3) {y - 2/4) Vv 2 - z ^v- 

27r Jy 3 
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Expanding the second square root for z\ small, each term in the series can be computed 
explicitly and it is most easily given in terms of a generating function [12] . 

F{a) = -TT^/(y 3 + a){y A + a) + |(y 3 + Vi + 2a) (3.24) 

as follows, 

oo 

Si = I (2/3 + 1/4)0/4 -?/3) 2 + ^X; c n A^W(0) 

n=l 

where c n are the coefficients in the expansion of yl — x and F^ n \a) is the n-th derivative 
with respect to a. 

The explicit answer has the following structure, 

S l = 9 -z 2 I-j- l K{z l) z 2) I 2 ), (3.25) 

where 

. 1/1 1 , , 2 1 \ 

K{x,y,z) = -xy 11 + — [x + y) + g^( x + V) + g^2 X?/ + ' " J ■ 

It is important to notice that this is symmetric in (x,y), namely, K(x,y,z) = K(y,x,z). 
This allows us to write, 

S 2 = -l Zl I+l-K{zi,z 2 ,I 2 ). (3.26) 

In the following we will simplify the expressions by putting m — g — 1. It will be useful 
to change variables to 

t = S 1 + S 2 , s = i(5!-5 2 ) (3.27) 
where t is the total 't Hooft parameter. Due to (13 .25j) and (I3.26|) one immediately obtains 

t = -{z 2 - z l ) y /l-2z 1 -2z 2 . (3.28) 

Note, that t can be regarded as a global parameter of the model. Different from t the ex- 
pression of s in terms of the Zi requires a transcendental function. This more complicated 
function reflects the dependence of s on the choice of the symplectic basis in (12. 2p . 

As mentioned earlier, there is another possibility to derive the periods as series in Z{ 
which was applied in [30]. There the authors consider a set of Picard-Fuchs differential 
operators, Ci,C 2 associated to the spectral curve and differential Q, which annihilate 
the periods. Therefore, these can be calculated as solutions to a system of ODEs. The 
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Picard-Fuchs operators, which are given in eq. flA.lj) of appendix |Aj have the following 
discriminant factors 

disc = z x z 2 I 2 J = z x z 2 {\ - 2(zi + z 2 ))(l - 6zi - Qz 2 + 9zf + Uz x z 2 + 9^). (3.29) 

Moreover, their solutions around z\ = and z 2 = describe the periods of the elliptic 
geometry (I3.18p . Due to the fact that one can find a combination of periods such that 
the mirror map becomes exact (I3.28p . it is convenient to introduce adapted coordinates 
Zi, i = 1,2, by 

Zi = Zi + z 2 , z 2 = -(zi - z 2 )\/l- 2(zi + z 2 ), (3.30) 
as well as coordinates ti, % — 1, 2, on the mirror by 

ti = s = i(5 1 -5 2 ), t 2 = t = 5 1 + 5 2 . (3.31) 

The Yukawa couplings may be found in eq. (1A.4|) as well as the genus one free energy 
F\ in eq. ( 1A.5I) . Due to the special type of the mirror map 

~z 2 = t 2 , (3.32) 

it is possible to derive a propagator which is of the following special form 

For the technical details as well as for the ambiguities that have to be computed we refer 
the reader to appendix |A] With the help of this input it is easy to implement the direct 
integration procedure for the cubic matrix model as outlined in section 12.11 It turns out 
that we can recursively construct the free energies up to genus four. Moreover, we can 
also evaluate F g (Si, S 2 ) for the cubic matrix model in perturbation theory, as was done 
in [34"| [30] . The expansions of our direct integration analysis read 

F2 = ~2io + if ) + f {Si ~ S2) + 338 ^ 2 _ 16325152 + 3385 2 + °( s3 ) 

Fs = iks (~h + ^l) + ^T^ Sl - S2) + f ( 525225 i 2 - 2734035iS 2 + 525225|) + 0(S 3 ) 

Fi = ~ lio (if + 5|) + ^T^ 1 - S ^ + IT ( 7883695 i ~ 43874365x52 + 7883695 2 2 ) + C(5 3 ). 

(3.34) 

These results agree with the low-order results obtained in [34"| 130] . In the following we 
explain how to parameterize the ambiguity and how to fix the unknowns entering our 
ansatz. 
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3.2 Direct integration, boundary conditions and integrability 

In the last section we set up the necessary ingredients to perform a direct integration of 
the holomorphic anomaly equations. As mentioned in section I2TT1 the free energies F g can 
be written in the following way 

3g-3 

F g =J2a k (z 1 ,z 2 )(S' z ^) k + f g (z 1 ,z 2 ), (3.35) 

k=l 

where a k are rational functions completely determined by the recursive procedure. f g 
is the holomorphic anomaly, which is not constrained by direct integration and must be 
fixed by supplying further boundary conditions. The amplitudes F g should be well-defined 
over the whole moduli space except for points at the boundary of moduli space where the 
elliptic geometry (13 . 1 8[) acquires a node, i.e. a cycle of S 1 -topology shrinks. Such points 
are known as conifold points and are given by the zero loci of the discriminant of the 
Picard-Fuchs system, which we also call conifold divisors. 

Thus, regularity and holomorphicity imply that f g should be a rational function of Zj, 
where the numerator is at most of the same degree as the denominator. The denominator 
is given by the discriminant factors and takes the form [z\Zi J 2 ) 29 ~ 2 . This gives the 
following ansatz for the holomorphic ambiguity 

(9) ^k A 



V n (9 'z k 
l^ki u k,i z i 



Z% 



where the a^?\ have to be determined by the boundary conditions. Note that due to the 
symmetry of the model in %\ and z 2 it is enough to restrict the numerator to a polynomial 
which is symmetric in Z\ and z 2 . In order to be well defined as z$ — > 00, the degree of 
this polynomial must be at most 12g — 12. It turns out that it is sufficient to truncate 
the degree at 9g — 9, as long as g < 4. However, this reduced ansatz may not be present 
at higher genus and one would have to deal with the full ansatz of degree 12(7 — 12. 

There are two boundary conditions which we will refer to in the following as a cycle gap 
and b cycle gap. Let us first consider the a cycle gap. Due to the Gaussian contribution to 
the partition function of a multi-cut matrix model (see for example [M] for more details) 
it is easily seen that the holomorphic expansion of F g at small filling fractions is of the 
form 

^2^Gp+5p) +o < s >- (9>1) - (M7) 
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Due to the absence of subleading singular terms in Si, this property of the expansion is 
referred to as the gap condition. The coefficients of the subleading singular powers of Si 
depend generically on the other Sj, with j ^ i -in fact they are (infinite) series in the Sj. 
Demanding the vanishing of these series leads in principle to an over-determined system, 
therefore in the multi-parameter case it is not easy to count the number of independent 
conditions implied by (13.37)) . 

The gap condition is also present in the expansion of genus g topological string am- 
plitudes near a conifold divisor [3J [291 12H], where we have 



Here, II is a flat coordinate normal to the divisor. In view of the Dijkgraaf-Vafa corre- 
spondence, this behavior should also characterize multi-cut matrix model amplitudes near 
the divisors of the spectral curve geometry. Again, since the coefficients of the subleading 
powers of II depend on the coordinates tangential to the conifold divisor, the counting of 
conditions in the multi-parameter case is not easily done. 

However, it turns out that, when both constraints, ( 13 .37)) and ( 13 .38)) . are taken into 
account, the holomorphic anomaly f g is completely and uniquely fixed. We checked 
this explicitly for genus g < 4. It is then natural to conjecture that the a and b cycle 
gap conditions are always sufficient to fix all unknowns in the holomorphic ambiguity 
for general matrix models with polynomial potential. Following [29] we refer to such a 
property as integrability of the holomorphic anomaly equation. 

3.3 Modular covariant formulation 

In the last sections we explained how to solve the cubic matrix model with the techniques 
known from topological string theory. However, we used a somewhat artificial description 
which does not make the symmetry properties of the geometry completely explicit. Such 
a formulation is given by writing all quantities in a covariant modular way. Since the 
geometry is an elliptic curve together with meromorphic differential Q we expect not only 
to parameterize the topological amplitudes F g by the elliptic modulus r but in addition 
by an auxiliary parameter. In the following we explore how this can be achieved in detail. 
We start by transforming the quartic curve ( 13 .18)) to Weierstrass form, where it is easy 





(3.38) 
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to read off the j-function. It is given by 

■ / v 16((l-3^-3z 2 ) 2 + 12^ 2 ) 3 

J(*l> **) = 77] 7 o A ^ ( 3 - 39 ) 

z\z 2 ((1 - 3^i - 6z 2 y - 4:Ziz 2 ) 
Comparing this modular invariant to its usual Fourier expansion 

j(r) = q" 1 + 744 + 196884g + 0(q 2 ), (3.40) 

we get a relation r = r(zi, z 2 ). Using the definition of j in terms of modular forms yields 
actually a rational expression. 

It is also easy to identify the auxiliary parameter which accompanies r. Note that the 
periods/filling fractions are taken with respect to the differential Q, which is meromorphic. 
Thus the sum of all cycles is of course homologically trivial, but the sum of the periods 
does not have to vanish and is rather proportional to the residue of Q. Since this residue is 
related to the auxiliary parameter, it is natural to parameterize the topological amplitudes 
by both, r and t = g s N. Due to (13. 28ft we obtain a relation t = t(z\, z 2 ). 

In principle this allows us to rewrite all quantities in terms of modular forms together 
with an auxiliary parameter t, by combining the rational expression r = t(zi,z 2 ) with 
t = t(zi, z 2 ). If we do so, we obtain 

u = 1 ~ 3Z 1Z. — , 4t= {z l -z 2 y\-2z l -2z 2l (3.41) 
where u is given in terms of modular forms b, c and d defined in appendix [C] by 

u = C -±±. (3.42) 
b 

However, it turns out that, for the general cubic matrix model, the resulting formulae 
become too complicated. The reason is that the corresponding spectral curve is a generic 
elliptic curve. However, if we specialize the calculation to the slice t — 0, or Si — —S 2 , 
the curve has T(2) monodromy (it is the Seiberg-Witten curve of [47]) and it is possible 
to exploit the formulation in terms of modular forms, as we will see in section HI 

Fortunately, it is possible to give some closed and simple expressions using modular 
forms for the genus zero and one sectors, which will prove to be useful in due course. We 
start by quoting the perturbative calculation of the planar free energy [12] 



Fo(Sl ,S 2) = l Sl log (iy + isl log (_^) - 3 (S , + sl) + 2SA log( m ) 

1/2 2 n (3.43) 

/ C3 cc2c 1 c C C2 c3 



^A 3 V3 1 1 2 3 2 



0(S 4 
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Of course, F (S\, S 2 ) is symmetric under the exchange S\ <H- —5*2. From the prepotential 
we can define the tau-coupling (12 .5ft and also introduce 



2vrir = (3.44) 



It was shown in [H] (see also [9]) that r can be computed in terms of elliptic functions as 

K' K(k') 

t =^='iM' (3 ' 45) 



where 



k = r dz = 2 tew, e - (X1 - X2)(X3 - X4) 



B1 y/\a{z)\ y/(x 1 -x 3 )(x 2 -x A ) Oi -x 3 )(x 2 -x A ) 

K! = I -^== = 2 = K(k') t k' 2 = 1 - k 2 . 



J-2 



y/W(z)\ y/{x\ ~ x 3 ) (x 2 - x 4 ) 



This modular parameter turns out to match with our definition by the j'-function men- 
tioned above. We find, in the full theory 

1 f-S 1 S 2 \ , 17(^-^2) 2(83^- 209^-1 + 83^) , 
7rir = - log — H 5 7 (3.47) 

Let us now consider genus one. Akemann [2] gave a simple expression for F 1 , that reads 
4 

F i = -24 E ln M * - 2 ln " 12 ln A + 8 ln(ar " ^ + 8 ln( ^ ~ at)2 ' (3 - 48) 
i=i 

where A denotes the discriminant of a(x). Using that Mi — g for the cubic matrix model 
as well as Thomae's formulae, cf. app. [Cj this can be written compactly as 

F x = -log V(r)-^ log A, (3.49) 
where 77 is the Dedekind eta-function. 

4 The cubic model for Si = —5*2 

In the following we will specialize the cubic matrix model studied in section [3] to the slice 
S\ = —S 2 . On this slice, t — 0, and the direct integration procedure simplifies. Moreover, 
we are able to write all quantities which are needed for direct integration in terms of simple 
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modular forms. The underlying reason for this is that, when t — 0, the spectral curve 
of the matrix model becomes the Seiberg-Witten curve, which has simple monodromy 
properties. Therefore the recursive procedure will be very efficient in obtaining results at 
high genus. 

First of all notice that, by contour deformation, 



, dz I 2m m 2 A a A 

Si + S 2 = g j> l + z + — z 2 + -z* + r^z* = — . 4.1 

2 4 W g g 2 g 2 g 2 2g 



2=0 

Therefore, if the parameter A in ( 13. 17ft vanishes A = 0, we have 

t = 5i + S 2 = 0. (4.2) 

In this case one also has [17] 

Tu = r 22 = -T12 = r, (4.3) 

which can be seen from (12.51) . 

From the point of view of the original matrix model, the slice Si = —S 2 involves an 
analytic continuation in the space of 't Hooft parameters. This is because on this slice 
Si/S 2 = Ni/N 2 = — 1 , which can not be implemented in the matrix integral ( 13.151) , since 
N lt2 are a priori positive integers. In terms of matrix integrals, the slice Si = —S 2 can 
be related to a super-matrix model [52l [201 EH]- A Hermitian supermatrix has the form 

(4 ' 4) 

where A (C) are Ni x Ni (N 2 x N 2 ) Hermitian, Grassmann even matrices, and ^ is a 
matrix of complex, Grassmann odd numbers. The supermatrix model is defined by the 
partition function 

Z s {Ni\N 2 ) = J V^e^ StlW ^\ (4.5) 

where we consider a polynomial potential W($) and Str denotes the supertrace. There 
are two types of supermatrix models with supergroup symmetry U(Ni\N 2 ): the ordinary 
supermatrix model, and the physical supermatrix model [52] • The ordinary supermatrix 
model is obtained by requiring A, C to be real Hermitian matrices, while the physical 
model is obtained by requiring that, after diagonalizing $ by a superunitary transforma- 
tion, the resulting eigenvalues are real. The partition function of the physical supermatrix 
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model reads, in terms of eigenvalues [521 [20] 

(4.6) 

where the two groups of eigenvalues fa, Vj are expanded around two different critical 
points of W(x). This partition function is related to ( 13. 15ft after changing N 2 — > —N2 
[20] . therefore it gives a physical realization of the S1/S2 < slice of the moduli space. 
Notice, that the moduli space of the local Calabi-Yau for generic complex Si, S 2 describes 
both the original matrix integral (I3.15P and its supergroup extension (14. 6p . 

4.1 The geometry 

In the following we discuss the geometry underlying the curve with A = 0. It is easy to 
see that, up to a shift in the x coordinate, it can be written as 

y 2 = (x 2 -a 2 )(x 2 -b 2 ), a>b. (4.7) 

If we compare this to the Seiberg-Witten curve (47] 

V 2 = (x 2 - uf - A 4 SW , (4.8) 

we find that they are equal once we identify the parameters as 

a 2 + b 2 , „ a 2 - b 2 . J . 

« = —3—, A iw = —j—- (4-9) 

We also want to translate these parameters in terms of the cubic matrix model variables. 
This was already done in [T7J [34] , and we have 

A = — , u = -A . 4.10 

We will set 

5 = 1- (4-11) 

On the other hand, we have the following relation between the A parameter appearing in 
( 13.431) and the Seiberg-Witten scale 

A = t^Asw- (4.12) 
20 



For the simple curve (14. 7p one can compute many quantities directly and relate them 
to modular forms or elliptic integrals. As a starting point the period integrals 



S — Si 



So 



U = d s F Q (S,-S) 



(4.13) 



can be computed in terms of simple elliptic functions, which was done for S in ref. |17j . 
Repeating this analysis yields 



s = ± 

2m 



y(x)dx = -J- (a 2 + b 2 )E(h) - 2b 2 K(k{ 

D7T L 



as well as 



n 



y(x)dx 



b 3 



[a 2 + b 2 )E(k[) + (b 2 -a 2 )K(k[) 



where the elliptic modulus k\ and its complementary one k[ are given by 



a 2 -b 2 



kl 



or cr 
The modulus k\ is related to the usual cross-ratio k 2 introduced in (I3.46P as 

Ik 



(1 + k) 2 

In order to obtain expansions of the periods we introduce the parameters 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



A 2 
u 



A 2 
u 



(4.18) 



Small /i corresponds to the semiclassical regime of Seiberg-Witten theory which occurs 
at it — > oo, whereas small //d relates to the region near u — > Ag W , where a magnetic 
monopole becomes massless. In these variables the periods read 



n 

3" 
M2 



S _ VTTJt 
4 V ^7I D " 



3 

«2 



E 



(l 



2/i 



E 



+ /i 



(a* - 1)^ 



(Md - 1)K 



2/i 
1 + Jl 



2 - /i D 



(4.19) 



(4.20) 



.2 - /i D . 

Note that S/u 3 / 2 and n/n 3 / 2 are dimensionless. Further, we expand ( I4.19P around /i = 
to obtain 

S fi 2 3/i 4 35/i 6 1155/x 8 45045/x 10 



3 

U2 



8 + 256 ' 8192 ' 524288 ' 33554432 



+ 



(4.21) 
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which is the expansion (4.19) of [33], after changing to the appropriate variables. The 
inverse expansion is given by 

2 / S \ 2 17/5 \ 3 375 / S \ 4 10689 / S \5 , jnn . 

" =8 7S-<^) "lb) - uKsaO -Tab) + "" (4 ' 22) 

We introduce now the following elliptic modulus r as 

r - i^M - i*^ - 1 ^ f4 23) 

which can be expanded in /i. By inverting this series one can derive /x as a function of To- 
In particular we observe 

/i=-^-, (4.24) 
c + a 

where we follow the notation^ of [30] . In turn the expression (I4.19P defines the variable fi 
as a function of 

8^ (4.25) 



u 3/2 m 3 

as well, and in particular the series (14. 23ft defines tq as a function of S/u 3 ^ 2 : 

, /S\ 5 B /S\ 2 71260 / S \3 
27rir = log — +34— + 750 — +——( — ) +--.. 4.26 

Moreover, comparing with (13 .47p yields the identity 

r = ^r(5,-5), (4.27) 

which is obvious also from (14.231) . 

Consider now the dual elliptic modulus t^d, obtained by a ^-transformation on the 
elliptic modulus, 

to,d = --• (4.28) 

To 

Following the same lines of thought as before, this defines t 0j d as a series in the dual 
period 

Sd " Sd (4.29) 



u 3/2 m 3 

In the following we will set 

Asw = 1 (4.30) 



For our conventions on modular forms used in this section, see appendix [C] 
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so in particular fi = u~ l = 4/m 2 . Note that (14.241) therefore defines m as a function of 
To. Strictly speaking, m is hence a function of S, but in order to establish the relation 
between tq and r(S,—S), i.e. (14.271) . we treated m as an independent variable. In all 
subsequent formulas and expansions we will do so as well. 
Next, we compute the Yukawa coupling 

C sss = (4.31) 
This follows from the general formula for two-cut matrix models given by [JT] 



ds 3 



Ah--- M 4 /C 3 Y[( Xi - x,) 2 ] ^[n^ n (x fc - x ; ) 2 ] (4.32) 



k<l 



where M, = M(xi), the spectral curve is written as in (]3.8p . and /C is given in (I3.46P 
When applied to the Seiberg-Witten curve (14.71) we obtain 

Csss = d ^ [To) = 64 v / 2( c + rf ) 5/2 - (4 _ 33) 
ds m 3 b 2 cd 

To see this, one has to apply Thomae's formula, which relates the branch points Xi to 
^-functions [2l] and further one has to express K, in terms of modular forms. This is done 
as follows. Note that 

f a dx 1 

K = / 77 2 Zf 2 » T = (434) 
Jb \f (a? - x 2 )(x 2 - b 2 ) a 



The dimensionless combination y/uJC can be expanded as a series in To since 
This yields 



u£ = -^^^ f ~~~ — ) . (4.35) 



We can check the formula f 14 . 3 1 D by calculating this quantity directly from the pertur- 
bative result. Evaluating the derivatives at Si = — S2 = S we obtain from (13.431) 

d 3 F 2 r m 3 5 S 2 1 . 

— ^ = i 34 + — + 1500— + 71260— + ■■■}. 4.37 

Using (14.261) this coincides with ( 14.31ft . if we treat m as an independent variable. 
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The expression (I4.33P for C sss is a modular form of weight —3 on the modular group 
r(2) defined in the Appendix. We will use as generators for the ring of modular forms on 
r(2), M*(r(2)), the functions 

K 2 = c + d, K 4 = b\ (4.38) 

which are modular forms of weight two and four, respectively. Note that instead of 
considering the T(2) description of the Seiberg-Witten curve (14. 7ft we could also use 
the equivalent T (4) description, which amounts to trade r for 2ro = r(S,—S) in all 
expressions of this section. 

4.2 Direct integration and higher genus amplitudes 

Having discussed the genus zero sector of the cubic matrix model specialized to the slice 
Si = — S2, let us now turn our attention to the higher genus free energies F g . According 
to [30j [22] the matrix model free energies F g can be promoted to modular invariant, non- 
holomorphic amplitudes F g (r , f ) which satisfy the holomorphic anomaly equations of [6] 
in the local limit. The matrix model F g is recovered by formally considering the limit 
fb — > 00 while keeping r fixed. 

In order to apply this, we must compute the full non-holomorphic genus one ampli- 
tude Fi and derive the propagator S ss . Using the general formula (I3.48[) specialized to 
the Seiberg-Witten curve ( 14. 7\\ . and by following the same argument as for the Yukawa 
coupling C sss , we obtain 

Fi{t , r ) = - log( v/l^WroM-fo)) + ~ log (^^) ■ ( 4 - 39 ) 
Indeed, when expanded we find 

Fl = _ 1 log S + S + 15 / S\ 2 + 6202 / S\ • + 32286 + . 

o dm 6 \m A J 9 ym- 3 J ym 4 J 

which is precisely the series for F\ obtained in [31] after setting S\ = —S2 = S. 
Next we turn to the propagator S ss , defined by 

Cl s = BsS ss , (4.41) 

where C S ss is the complex conjugate of the Yukawa coupling C sss and the indices are 
raised by means of the metric 

G sS ~ Imr . (4.42) 
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Using the chain rule and the relation (I4.3ip yields 



^(to, r ) = -^C sss E 2 (t , f ) + dji(r ), (4.43) 



where fi is given by 

fi(Tn) = 

4 \VK, 
Hence, the propagator is identified with 

S ss = -±E 2 (T ,f ). (4.45) 

Now we are prepared to apply the method of directly integrating the holomorphic 
anomaly equations according to [26], |3]. In the conventions of this section the holomorphic 
anomaly equations can be cast into the following form 



dE 2 ~ 192 ° ss 



D T C sss - 9-1 
Dt F 9-~l 7} sss D T0 F g ^i + y y j D T0 F h D n) F g -. h 

^sss , , 

h=l 



, (<?>1) (4.46) 



where D T0 denotes the Maass derivative acting on (almost-holomorphic) modular forms 
of weight k as 

D T0 = — — - . (4.47) 

27ri dro 47rImro 

Since the ring M*(r(2)) = C[E 2 , K 2 , K4] is closed under D To , and the F 9 's are modular 
invariant forms, the holomorphic anomaly equation can be integrated with respect to E 2 . 
We obtain the following schematic result 

3g-3 

F 9 (r , 7b) = A 2 " 29 • £ c[ 9 \r )E k 2 (T , f ) + f g (r ), (4.48) 
fe=i 

where c^(tq) are modular forms of weight 8(g — 1) — 2A;, completely determined by the 

holomorphic anomaly equation, and A is just the denominator of C sss . In particular it is 
a weight eight form given by 

A = m z (Kl - K A )K A . (4.49) 

All the non-trivial information is encoded in the holomorphic ambiguity / 3 (to). It has to 
be derived genus by genus by supplying further boundary conditions. In the particular 
case of the cubic matrix model specialized to the slice Si = —S 2 , we will argue in the 
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next subsection that / s (to) can be fixed at all genera. Applying this procedure we were 
able to integrate the holomorphic anomaly equations and obtained the matrix model free 
energies to genus 52. 

Let us at least present the result for the full non-holomorphic genus two amplitude 



160Af 93 16K*(5K 2 - 7K 4 ) 



-2 



F2(To ' To) 81m 6 (Af - KtfK'f 2 9m e (K 2 - K^K 2 ^ 



1 



8Af (77A 2 4 - 132K 2 K 4 + 63iff) ~ 4A|(2051A 2 4 - 4005Af A 4 + 1890ATJ 



27m 6 (Af- A 4 ) 2 Af z 405m 6 (Af - A 4 ) 2 Af 

(4.50) 

Here we collect some low genus expansions of the free energy amplitudes of the cubic 
matrix model on the slice Si = —S 2 = S: 

„ 1 m 6 35 S S 2 1341064 S 3 S 4 

m = ~u 0l p + y ~, + ^08-, + -^--y + 24734074-^ + • • • 

12 1 m 12 10010 5 4036768 S 2 1883381692 S 3 S 4 

m 12 F 3 = —-^r + — 5 + ^ F + z q + 38608040638—^ + ■ ■ ■ 

504 S 4 3 to 3 3 m 6 7 to 9 m 12 

18 r. 1 ml§ 8083075 S ^ AnAMnAn S 2 4618613451580 S 3 

m 18 ^4 = -zttt^t + ~ + 1749491040— 7 + 5 + • • • (4.51) 

720 S 6 3 m 3 m 6 9 m 9 

94 1 m 24 1 3013750750 5 S* 2 17515677810823140 S 3 
™ F = = 558^ + 3 J? + "^8280413440^ + jj ^ + ■ ■ - 

m » F . _ _ J" + 1169936192 42504 + «""20800 S1 + ■ 
163800 o iu m 3 3 m° 

We can check some of these results by comparing to the perturbative calculations of 
[34] specialized to Si = —S 2 = S. We observe agreement for genus two and three at low 
order in S/m 3 . All higher genus computations are new results. 

The direct integration procedure outlined here is by far the most efficient method to 
calculate higher genus amplitudes in matrix models. It only takes a few minutes to reach 
e.g. genus 10 on a conventional personal computer. 

4.3 Boundary conditions and integrability 

According to [U [22] F g is an almost-holomorphic modular invariant form under the space- 
time duality group, in this case T(2). Hence, F g is regular except for some points on the 
boundary of moduli space. 

Regularity and holomorphicity imply that f g should be a rational function, where its 
denominator is given by an appropriate power of the discriminant of the curve. From the 
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expression (14.481) we see that the denominator of f g is given by A 29 2 , hence a weight 



8(g — 1) form. Modularity now implies that the numerator has to be a form of finite 
weight, in order to cancel the weight from the denominator. Since the space of weight 
k forms is finite dimensional, there are only finitely many coefficients to determine. In 
particular, for T(2) we have 



In summary this justifies the ansatz 

4(9-1) 

f g (r ) = A 2 ' 29 ■ a k KfKf 9 - X) -\ (g > 1) (4.53) 

where A is given in eq. (I4.49p . This implies that there are 4g — 3 unknown constants a k 
in the ambiguity f g . These are completely and uniquely fixed by imposing the following 
two boundary conditions. 

First, we know that the holomorphic expansion of F g at small S has the structure 
(13.371) specialized to the slice, which imposes 2g — 1 conditions on f g and leaving 2g — 2 
unknowns. Further the holomorphic expansion at conifold divisors is of the form (I3.38|) . 
where II is a suitable coordinate transverse to the divisor which vanishes at the conifold. 
In our case II is the dual period. Thus, (13.381) imposes 2g — 2 further constraints on the 
ambiguity, and it determines it completely. 

5 Non-perturbative aspects 

In this section we address non-perturbative effects in the two-cut matrix model, and its 
connection to the large order behavior of the 1/N expansion. We first review the one-cut 
case. 

5.1 Non-perturbative effects in the one-cut matrix model 

For concreteness, we will focus here on the cubic matrix model which we are analyzing in 
this paper. 

In the one-cut cubic matrix model, the large N limit is described by a distribution 
of eigenvalues around the minimum of the potential at x — 0. The eigenvalues fill the 




k > 2, 
else. 



k even. 



(4.52) 
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interval [a, b]. It has been known for some time that instanton sectors in this model are 
obtained by tunneling a finite, small number of eigenvalues £ <C N from this interval 
to the maximum of the effective potential, located at xq. The structure of the partition 
function in the £- instanton sector has been determined in HI] , and at one loop it has 
the form 2 

Z {£) = -^ m G*(£+l)rf exp f-^j |l + C%)}. (5.1) 

In this equation, G2(z) is the Barnes function. A is the instanton action, and it can be 
computed in terms of the spectral curve of the one-cut matrix model as 

A= dzy(z). (5.2) 
Jb 

Finally, \x\ is the one-loop contribution, and it has the explicit expression 

b-a 1 

fix = — (5.3) 

^M(x )[(a-x )(b-x p 

In [ID] it was argued, following standard arguments in the large order behavior of per- 
turbation theory [27], that the free energy of the one-instanton amplitude, should 
determine the leading asymptotics at large g of the perturbative amplitudes F g , according 
to the formula 



If we write 



oo 



F^=9l /2 e- A ^J2^t\ (5.5) 
we obtain the full 1/g asymptotics 

1 A-2a-bT^(n„ , i\r UlA 1 ' 1 



F 9 - 9 -A~^T(2g + b) V 



where 

5 



(5.6) 



(5.7) 



2 

The formula (15. 6p can be regarded as a generalization of the asymptotics for formal 
solutions of nonlinear ODEs. The reason is as follows. In the double-scaling limit of the 
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matrix model (see [IS]), the total free energy of the matrix model becomes a function of 
a double-scaled variable z, 

F(t,g s )^F ds (z), (5.8) 
and the specific heat u = —F^ s (z) satisfies the Painleve I equation 

u 2 - ^u" = z. (5.9) 

In particular, the genus expansion of the cubic matrix model leads to a formal solution 
of Painleve I 

oo 

u(z) = z 1/2 J2 u 9fi z ~ 59/2 - (5-10) 

9=0 

On the other hand, the instanton sectors of the matrix model lead to instanton corrections 
of the form 

oo 

Ul {z) = z V2-M/8 e - iaz5/4 u n/ z- 5n/4 (5.11) 

n=0 

where 

a = 8 -^. (5.12) 
5 

It can be shown that the coefficients of (I5.10P have an asymptotic behavior at large g 
which is governed by the one-instanton solution U\(z) in (15.111) . The precise formula is, 

a~ 29+ Kf~ ^ Una i 











7T1 [ 



«*o ~ 9 — r2 S - ~ 5 l + V . „ . (5-13) 



2J 



where S\ is a Stokes constant. One can explicitly check [TUHO] that ( 15 . 1 3[) can be deduced 
from the double- scaling limit of the asymptotics (15. 6p . In particular, the constant a is the 
double-scaling limit of the instanton action. 

5.2 Non-perturbative effects in the cubic matrix model 

Non-perturbative effects in multi-cut matrix models have been studied in [91 ST]. A 
multi-cut matrix model with a fixed choice of filling fractions must be regarded as a fixed 
background, and any other choice of filling fractions leads to an instanton correction to 
the free energy on the fixed background. To be concrete, let us consider a two-cut matrix 
model with a fixed background given by N\, N 2 eigenvalues in the stable and unstable 
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saddle points, respectively. The partial 't Hooft parameters Si, S2 are given as usual by 
Si = g s Ni. The total partition function is of the form 

Z = Z(N U N 2 ) + ^ C' Z ( N * -*,N 2 + £). (5.14) 

The sum over i corresponds to the tunneling of i eigenvalues from the first cut to the 
second cut, and at large N, the corresponding partition functions have the form 

z (l) = C y 2 /2 exp (^JA^j |i + (g s )y £ G IT (5.15) 

where 

A = d s F and g = exp(d s 2 F ) . (5.16) 

The variable s is given in (13.271) . If the cuts of the matrix model are the intervals [x±, X2], 
[23, £4], the instanton action A can be written as 

A= y(x)dx. (5.17) 

If Ke(A) 7^ 0, the instanton contributions are exponentially suppressed if sgn(Re(A)£) > 0, 
and they are exponentially enhanced if sgn(Re(v4)£) < 0. This is just reflecting the fact 
that the generic background is unstable and if we expand around it we will find tachyonic 
directions. Notice however that both corrections are non-perturbative in g s , therefore 
they are invisible in the genus expansion. 

It is generically expected that the existence of these non-perturbative sectors leads to 
the factorial divergence of the genus expansion around a fixed background. The growth 
of the perturbative string amplitudes at large genus (and fixed Si, S2) should be of the 
same form as in (I5.6p . i.e. 

i^(Si, S 2 ) ~ 9 A~ 2 °- b T(2g + b) + Oig- 1 ) (5.18) 

where A is given by ( I5.17P and b is a constant. 

We can test these predictions by numerical methods using our results from direct 
integration. We start by concentrating on the slice Si = — S2. Note that in this case 
the instanton action A is given by the dual period II, whose explicit expression is given 
in eq. f )4.15p . In order to extract the asymptotic of the sequence {Fg} g > we employ a 
standard numerical technique known as Richardson extrapolation. The method removes 
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the first terms of the subleading tail and hence accelerates the convergence. Given a 
sequence {S g } g >o in the form 

S fl = a + - + ^ + ... (5.19) 

9 g 2 

its Richardson transform is defined by 

M 9 ,N) = E i ~ 1 tZ { -tv )N s ' +k ' (5 - 20) 

fc>0 ' v '' 

such that the sub-leading terms in {S g } g >o are cancelled up to order g~ N . In fact, it can 
be shown that if {S g } g >o is a finite sequence, the Richardson transform returns exactly 
the leading term do- 

Comparing (I5.19P with ( 15 . 1 8[) one can extract the instanton action by considering the 
sequence 

F s+l ^ 1 / 1 + 26 _ 2 , 



Q ^wFr*v + — +oi9 } )- (521) 

Once A is confirmed, one can then obtain the parameter b from the new sequence 

Q' g = 2g (A 2 Q g - 1) = 1 + 2b + 0{g~ l ). (5.22) 

In Fig. [2] and Fig. [3] we plot the sequences Q g , Q' g , together with their Richardson trans- 
forms, for two values of S. It is obvious from the numerical calculation that the large 
genus asymptotics is controlled at leading order by the instanton action. In addition, we 
find numerically that 

b = -l. (5.23) 

This value of b is different from the one characterizing the one-cut model (|5.7p . In fact, 
the value ( 15.71) corresponds to the universality class of pure two-dimensional gravity, while 
the value (15.231) corresponds rather to the universality class of the c = 1 string [45]. It 
is interesting to see that both behaviors are present in the two-cut cubic matrix model, 
along different submanifolds of the moduli space (the 2d gravity behavior takes place in 
the slice S2 = 0, while the c = 1 behavior takes place in the slice Si + S2 = 0) . 

Turning our attention to a generic value of the filling fractions (Si, S2) in the cubic 
matrix model, we can try to test our prediction (I5.18P by using the results from direct 
integration of section [3] Since we computed F g up to genus four, we can only explore the 
first four elements of the sequence {Q g } g >o, eq. ( 15.211) . In order to have a better control 
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Figure 2: The sequence Q g (•) and two Richardson transforms (■, ♦) at To = \ (left) and 
(right) which corresponds to S ~ 0.139 and S ~ 0.117 + 0.016i, respectively. 



r 



3+9 



The leading asymptotics as predicted by the instanton action \A\ is shown as a straight 
line. The error for genus 52 is about 10 -8 % and 10 -10 %, resp. 











Figure 3: The sequence Q' g (•) and two Richardson transforms (■, ♦) at To = ^ (left) and 
r = j + | (right) which corresponds to S ~ 0.139 and S* ~ 0.117 + 0.016i, respectively. 
The leading asymptotics as predicted by the parameter b = —1 is shown as a straight 
line. The error for genus 52 is about 10~ 8 % in both cases. 



of the error, we consider a perturbation around the submanifold Si + S2 = where the 
large order behavior is well established. The instanton action ( I5.17P is calculated using 
(OJ) by 

dF dF 



_ 

dS! dS 2 
= log(Si)Si - \og(S 2 )S 2 - 



(5.24) 



S 2 + 0{S 7 



Fig. H] shows Q3, Rq(1,2) and \A\ 2 as a function of Si in the vicinity of the slice point 
Si = — S 2 = S = 0.004, where convergence is ensured. We observe that the behavior of 
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Qs and Rq(1, 2) is qualitatively the same as predicted by the instanton action. Moreover, 
their relative errors stay roughly constant over the complete data set. This seems to 
indicate that the large order behavior of the genus expansion is also governed by the 
instanton action in the general two-cut cubic matrix model. 

Unfortunately, our numerical results for the generic case are not good enough to deter- 
mine the value of b reliably. It is an interesting question to know how this value changes 
as we move in the moduli space. We expect it to be b = — 1 except in the one-cut slices 
Si = 0, S 2 — 0, where it takes the value (15.71) . 




Figure 4: Q 3 (•), Rq (1,2) (■) and \A\~ 2 (♦) are plotted for several values of Si around 
the slice point Si = —S 2 = S = 0.004. Q 3 and Rq(1, 2) have a relative error of about 30 
% and 10 %, respectively, as compared to the instanton action \ A\~ 2 throughout the data 
set. 



5.3 Asymptotics and non-perturbative sectors 

In principle, one should be able to refine the asymptotic formula ( 15. 18ft and obtain a 
generalization of (I5.6P involving the g s expansion of instanton solutions. A natural guess 
is that the relevant instanton solutions are the closest ones to the given background, i.e. 
the instanton amplitudes ( I5.15P with I = ±1. This guess would relate the large genus 
behavior of F g (ti,t2) to an integral of the form ( 15. 4ft . involving this time and F^ 1 ^. 



33 



However, this expectation turns out to be too naive. Indeed, it seems that the asymptotics 
involves new non-perturbative sectors whose matrix model interpretation is yet unknown. 



Figure 5: The asymptotics of the coefficients of the £-th instanton solution ut(z) of 
Painleve I is determined by the two nearest neighbor instantons, which are obtained 
by eigenvalue tunneling, and by the generalized instanton amplitude u^i, which is repre- 
sented here by the label (£,!)■ 

In order to explain this in some detail, we will come back to a simplest case where the 
asymptotics can be fully determined, namely the Painleve I equation and its instanton 
solutions ue(z). It is natural to ask what is the asymptotics of the coefficients u n ^ appear- 
ing in (15. lip . Notice that, when i is big, this instanton solution is the double-scaled limit 
of a two-cut solution, therefore the question of the asymptotics of this sequence is closely 
related to the original question concerning the asymptotics (15.181) . We have seen in (15.131) 
that the asymptotics of the perturbative solution is governed by the one-instanton solu- 
tion. In the same way, one would think that the asymptotics of the ^-instanton solution 
is governed by the £ ± 1 instanton amplitudes. It has been shown in [25] that this is not 
the case. In order to understand the asymptotics of a generic instanton sector, one has 
to consider more general amplitudes, labelled by two non-negative integers: 



The amplitude where m = is the standard instanton amplitude: ue\o{z) = ue(z). The 






(5.25) 
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other amplitudes can be obtained by requiring 

u(z,C h C 2 ) = Yl u n]m (z)C?C? (5.26) 

n,m>0 

to be a formal solution to the Painleve I equation, for arbitrary C±, C 2 , and that 

u n \ m (z) ~ e -("- m ^ 5/4 , z^oo. (5.27) 

The two-parameter solution of the Painleve I equation (15.261) is called a trans-series so- 
lution, and it was introduced by Jean Ecalle in the context of resurgent analysis (see 
for example [16] for an simple introduction to resurgence). It turns out that the asymp- 
totic behavior of the coefficients u n ,i in the £-th instanton ug is governed by the solutions 
U£±i(z), but also by the solution Um(z). This means that the asymptotics of the coeffi- 
cients u n> £ as n — > 00 can be obtained by a relation similar to (15.41) . but involving ue±i(z) 
as well as Ufu{z). For example, let us consider the one-instanton solution, and let us ask 
what is the asymptotics of the coefficients u ns i appearing in ( 15.111) with i = 1. An analysis 
based on resurgence theory, which can be verified with Riemann-Hilbert techniques, leads 
to the formula [25] 

«„,: ~. »-" +1,2 §-r(n - l/2){2«,, + (-1)>„, 2 + ± g^i±i=lg^K } (5.28) 

where u n ^ are the coefficients of the two-instanton expansion (£ = 2) in ( 15. lip , and /i n 2 
are the coefficients of the function 

«i|iW=^ /4 ^W 5n/4 (5.29) 

n>0 

It can be seen, by plugging (I5.26P in the Painleve I equation, that this function satisfies 
the linear inhomogeneous ODE 

- q u i\i + 2«o«i|i + 2ti 1 -u |i = 0. (5.30) 

There are similar, but more complicated, formulae for the asymptotic behavior of the 
coefficients u n / for arbitrary £, see [25]. They all involve the trans-series solutions ue\i(z). 

The instanton amplitude u^ +1 (z) can be obtained from the solution u^z) by tunneling 
one extra eigenvalue to the unstable saddle, while the amplitude u^i{z) can be obtained 
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\ / 




Figure 6: In the two-cut case, given a background with perturbative amplitudes F g (t 1 ,t 2 ), 
there are two instantons which are obtained by eigenvalue tunneling, and two generalized 
instanton amplitudes represented by (JVi, 1), (N 2 , 1). 

from ug(z) by tunneling one eigenvalue back to the stable saddle. The amplitude u»i(z) 
does not seem to have, however, an eigenvalue interpretation of this type. The different 
non-perturbative sectors governing the asymptotics of the £-th instanton solution are 
depicted in Fig. [51 

We can now come back to the original problem of determining the large order behavior 
of F g (Si,S 2 ), corresponding to a two-cut model with Ni, N 2 eigenvalues around the 
two saddles. There are two instanton configurations which are obtained by eigenvalue 
tunneling, with fillings N\ ± 1, N 2 =F 1- It can be seen that the sub leading terms in 
the asymptotics are not reproduced by using just these two configurations. This is not 
surprising, in view of the result for the instantons of Painleve I. The above analysis, based 
on [25], suggests in fact that there should be two other non-perturbative configurations 
in the two-cut matrix model, that we denote by (N\, 1) and (N 2 ,l), in analogy with 
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our notation in Fig. |5j These configurations are depicted at the bottom of Fig. [6j The 
instantons obtained by eigenvalue tunneling can be easily calculated from the free energy 
of the generic two-cut matrix model. However, we do not know how to compute the 
amplitudes involving these new configurations, since there is no analogue of the Painleve 
I equation (or the pre-string equation) for the generic two-cut matrix model. In general, 
it seems that the most general saddle-point of the two-cut matrix model should be labeled 
by two pairs of integers, (A r i,Mi), (A^ 2 ,M 2 ), associated to the two critical points of the 
cubic potential. 

6 Conclusions and further directions 

In this paper, building on [30], we have shown that the direct integration of the holo- 
morphic anomaly equations provides a powerful tool to calculate the 1/N expansion of 
multi-cut matrix models. We have seen that, in some circumstances, we can easily fix 
the holomorphic ambiguity and obtain explicit expressions for the genus g amplitudes. 
In general, we expect the anomaly equation to be integrable, in the sense that the gap 
conditions completely fix the holomorphic ambiguity. In the case of the two-cut cubic 
matrix model in the slice Si = —S2, we can use this method to determine the amplitudes 
to very high genus. 

These high genus results have allowed us to obtain quantitative evidence for the con- 
nection between large order behavior and eigenvalue tunneling in a multi-cut matrix 
model. However, our results indicate that the detailed large genus asymptotics of the am- 
plitudes cannot be understood just by considering the non-perturbative sectors associated 
with eigenvalue tunneling. Indeed, in a similar asymptotic problem analyzed in [25], it 
was necessary to include new non-perturbative sectors. It is only natural to suggest that 
a correct understanding of the asymptotic properties, in the multi-cut case, requires also 
the inclusion of new non-perturbative sectors. In the one-cut case and its double-scaling 
limit, the amplitudes in these new sectors can be obtained algebraically, as trans-series 
solutions to the pre-string equation and the Painleve I equation, respectively. In the 
multi-cut case there is no analogue of these equations, and therefore the corresponding 
generalized amplitudes can not be computed with our present tools. 

One obvious question is then the following: what is the interpretation of these new 
non-perturbative sectors in terms of matrix models or topological strings? We will give 
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now some hints which might help in answering this question. Let us first discuss the trans- 
series solutions u n \ m (z) appearing in (15. 26ft . It turns out that Uq\i(z) can be obtained from 
U£\o(z), the standard instanton amplitude, by changing the sign 

This corresponds to changing the sign of the string coupling constant g s — > —g s . If we 
think about the ue(z) as describing a set of i D-branes, then the natural interpretation of 
Uq\i (z) is as a set of I anti- D-branes. Indeed, it has been argued that anti- D-branes are 
obtained from D-branes in topological string theory just by changing the sign of the string 
coupling constant [50]. More generally, these should be the ghost D-branes introduced in 
p[3] . which reduce to anti-D-branes in the topological string context. It is then natural to 
interpret the generalized instanton amplitude u n \ m (z) as representing a state of n D-branes 
and m anti-D-branes at the unstable saddle, in the background of N — n + m D-branes 
in the stable saddle. If this interpretation is correct, the generalized amplitudes in the 
multi-cut matrix model, which we labeled by two pairs of integers (N\,Mi), {N2,M-i), 
should correspond to a saddle where there are Ni branes and Mi anti-D-branes at the z-th 
critical point, i = 1,2. 

One problem with this interpretation is that, as argued in [501 [19], such a configuration 
is described in principle by a quiver or supergroup matrix model. If this is the case, the 
non-perturbative configuration characterized by (N i: Mi), i — 1,2, would be equivalent to 
a configuration with only branes or only antibranes at the critical points. More precisely, 
we would get \Ni — Mi\ branes or |jVj — Mj| anti-branes depending on the sign of Ni — Mi. 
Since explicit calculations show that the amplitude u n \ m {z) is not equal to the amplitude 
u n -m\o{z) [25], the interpretation in terms of brane/anti-brane systems might not be 
completely appropriate. 

We believe that the appearance of these new sectors indicates that we do not fully 
understand the non-perturbative structure of matrix models and of two-dimensional grav- 
ity. Therefore, it would be very important to clarify their meaning and to compute their 
amplitudes in the multi-cut case. 
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A Data of the two-cut example 

In the following we collect the necessary data for our two-cut cubic model of the main 
body text. We restrict ourselves to the points in moduli space which are relevant for 
our discussion. For further background on e.g. the monodromy around several divisors in 
moduli space we refer the reader to [30] . 

A.l Large Radius 

d n c 2 = { Zl = 0} n {z 2 = o}: 

The Picard-Fuchs operators governing the periods of the cubic matrix model are given by 

d =(3 - 2zi - 6^ 2 )<9i - 2zi(l - 2z x - 6z 2 )d% + (1 - 10*i + 12*? + ^z lZ2 )d x d 2 
+ (3 - 6zi - 2z 2 )d 2 + (1 - 10z 2 + Az x z 2 + Ylz\)d x d 2 - 2z 2 (l - 6z t - 2z 2 )d\, 

C 2 = - 3(1 - 12zi + 18zf + 14^ 2 ) + {-3z 2 (l - 3z 2 + 2z\) 
+ zi(7 + 46^ - 18z 2 + 2<qz\ + zi(-39 + §2z 2 )))d x 
+ (-1 + 2z x + 2z 2 )(-2z x (l + 52? - 2z x z 2 - 3z\ - A( Zl + z 2 ))df 

(A.l) 

+ (zi + z 2 )(l - 8 Zl + 6zl - 6z 1 z 2 )d 1 d 2 ) 

- 3(1 - 122 2 + 14zi^ 2 + 18zl) + (-3^1 (1 - 3zi + 2z\) 
+ z 2 (7 - 18 2l + 26zl + (-39 + 62 Zl )z 2 + 46z 2 2 ))<9 2 

+ (-1 + 2z x + 2z 2 )((;z 1 + z 2 )(l - 8^ 2 - G Zl Z2 + 6z 2 2 )<9 1( 9 2 

- 2z 2 (l - 3z\ - 2z x z 2 + 5z 2 2 - A( Zl + z 2 ))d 2 2 ). 

Its discriminant can be determined to be 

disc = z x z 2 l 2 J = z x z 2 {\ - 2( Zl + z 2 ))(l - 6 Zl - 6z 2 + 9z\ + Uz t z 2 + 9zj), (A.2) 
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and its solutions around = 0, % — 1, 2, are given by the following expansions 

£i 1 

51 = — - -zi(2zi + 3z 2 ) + ... 

5 2 = -^ + -z 2 {3 Zl + 2z 2 ) + ... 

4 o 



(A.3) 



(A.4) 



IT = S 1 log + — (2z? - Wz lZ2 - hz 2 2 ) + ... 

n 2 = S 2 log (-|) - 1 + | - 1(5*? + 10^z 2 - 2 2 2 ) + . . . . 

The Yukawa couplings are given by 

1 - 6zi + 9z 2 - 5^ + 9^1 z 2 + 6z| 
212121 ~ 16z^P 

1 - 32j - 5z 2 

l-5z 1 - 3z 2 

Lj z 1 z 2 z 2 — -^gj 2 

1 - 5^i + 6z 2 - 6z 2 + 9ziz 2 + 9z% 
222222 ~ lQz^P ' 

where all other combinations follow by symmetry. The genus one free energy can be 
written as 

Fi = ~ log (det(Gy ) - 1 log(zi^) - ~ log J + i log J. (A.5) 
It is convenient to introduce new variables Zi, i = 1,2, by 

zi=zi + z 2 , z 2 = -(zi - z 2 )^l- 2(zi + z 2 ), (A.6) 

as well as coordinates tj, i — 1, 2, on the mirror by 

ti = s = ^(Si-S 2 ), t~ 2 = t = Si + S 2 , (A.7) 

such that the mirror map becomes as simple as possible. E.g. we have that 

z 2 = t 2 - (A. 8) 

This implies that some of the Christoffel symbols vanish: 

rgs, = «, for i = l,2. (A.9) 
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There are only four non- vanishing ambiguities of equation (I2.12p . that are given by 



1 i.i 

K - 98?, 4- K9.5 



5 - 285i + 52zf - 32zf - 112% 
2(1 - 2z 1 ){l - 85i + 20S 2 - Wzf + 16z 2 



~2\ ' 



pi 24 ^ 2 fA in) 

1/5152 1 - 85i + 205 2 - 16zf + I62 2 ' 1 ' 

Pl _ s - m, 

J z- 



2222 1 - 8zi + 20Sf - 165? + 165| ' 
This results in a propagator that has one non-vanishing component in 5-coordinates, i.e. 

3*1*1 = 45 2 _ g 45 2 + 44 ~3 _ 8325l5 2 + _ _ _ ? 5*1*2 = S 'z 2 Z! = S 'z 2 z 2 = Q ( A n ^ 

The covariant derivative closes on this propagator when one fixes yet another ambiguity 
, cf. eq. (12.111) . The only relevant, non- vanishing component is given by 

2l2l = 8(1 - 25 1 ) 3 (5 1 - A~z\ + 45? - 645j + 1445x2*) 
/zi (1 - 821 + 205 2 - 165? + 165 2 ) 3 ' 1 ' ' 

A. 2 Conifold 

Conifold J = {1 - 6z 1 - 6z 2 + 92 2 + 142i2 2 + 92 2 = 0}: 

We consider the point (zi,z 2 ) = (|, |) G J. Convenient coordinates are given by 

z c ,i = — t=(zi - z 2 ), z Ct2 = 1 - 4(zi + z 2 )- (A. 13) 

a/2 

z Ci i parametrizes the tangential direction to the conifold divisor, whereas z Ct2 the normal 
one. Transforming the Picard-Fuchs system to these new coordinates the polynomial 
solutions are given by 

1 



Ui — z c,iy 1 + z c,2 — z c ,i + 2^,1^,2 - g z c,iz c2 + 0(z c ), ^ 
u 2 = zlz + Szl^ + Oizi). 
We choose as flat coordinates 

t c>i =Ui, i = l,2. (A.15) 

By Inverting the above relations it is easy to calculate the holomorphic limit of the metric 
and the Christoffel symbols in z c coordinates. Transforming the Yukawa couplings 
as well as the ambiguities yields the propagator at the conifold point. This allows now 
to expand the free energies F g in the holomorphic limit at the conifold point. 
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B The one-cut solution 



A special case of the multi-cut matrix model occurs when all the partial 't Hooft param- 
eters are zero except for one, S2 — ■ ■ ■ — S n — 0. This is called the one-cut matrix model. 
A powerful, recursive solution of the one-cut matrix model at all genera has been known 
for a long time [7J 18] , and it is based on the technique of orthogonal polynomials. In this 
appendix we list some ingredients of the one-cut solution. 

We first recall that the orthogonal polynomials p n (A) for the potential V(X) are defined 

by 

/dA 
— e as p n (\)p m (X) = h n 5 nm , n > 0, (B.l) 

where p n are normalized by requiring that p n ~ A n + • • • . It is well known (see for example 
|16j ) that the partition function of the matrix model can be expressed as 

N-l N 

Z N =l[h l = h»l[r?~\ (B.2) 

i=0 i=l 

The coefficients 

Tn = ~7~~~ (B.3) 

satisfy recursion relations depending on the shape of the potential. They also obviously 
satisfy 

r n = (B.4) 

In the limit N — > 00, n/N becomes a continuous variable that we will denote by z, and r n 
is promoted to a function, R(z,g s ). The perturbative g s expansion is obtained by writing 

00 

R(z,g,) = ^2Re(z)g^. (B.5) 
We now consider the cubic matrix model with potential 

1 .W(z) = — (-z+^]. (B.6) 



9s 9s V 3 

The resulting matrix model is equivalent to the one we considered in the bulk of the 
paper, up to a linear transformation z — > az + b, but leads to a simple recursion relation 
for the coefficients r n (see for example |16j ) 



(yi-r n -r n+1 + a/1 - r n - r~Tj = g s n. (B.7) 
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The continuum limit of equation (IB.7P above is 

R(t, g s ) (y/l-R{t,g a )-R{t + g a ,g s ) + y/l - R{t,g 3 ) - R{t - g 8 ,g 8 j) = t. (B.8) 
At lowest order in £ and g s we find 

2R (t)^l-2R (t)=t. (B.9) 
It turns out to be convenient to express everything in terms of 



U(t,g s ) = s/\-2R{t,g a ) (B.10) 
This has a g s expansion of the form 

U(t,g s ) = J2m)9s 9 - (B.ll) 

Notice that 

u = U (t) (B.12) 

satisfies 

u(l-u 2 )=t, (B.13) 
and one chooses the solution with expansion 

t S/ 2 

« = i-5-T-- (R14) 

We also have the inverse relation 

R(t,g s ) = ^(l-U 2 (t,g s )). (B.15) 
This variable eliminates square roots. The recursion becomes 



U(l -U 2 ){Jl + (^-Y + Jl+ (^-) 2 \= 2V2t (B.16) 



u ) v \u 

where 

U ± = U(t±g s ,g s ). (B.17) 

Expanding the equation (IB. 16j> in power series of g s , and solving it recursively, one can 
compute in this way the coefficients U n (u), and from them derive the coefficients Rg. The 
first few read, when expressed in terms of 

1 - u 2 

r = - JL (B.18) 
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as 



Ri(t) 



(9r - 5) 



R 2 {t) 



3 (162r 3 + 1017r 2 - 1316r + 385) 



32(1 -3r) 4 ' 2048(3r-l) 9 
Once R(z,g s ) are calculated, one can easily calculate the total free energy 

F(t,g s ) = Y,9 2 s 9 - 2 F g (t). 



Let us define 



One then obtains [7J |8]: 



[z,g a 



R(z,g s 



ft R rl 2 ?" 1 

9 ]F = / d z (t - Z ) log S (z) + £ 9? tS, £^ 







(t - z) logZ(z,g s 



z=0 



+ 



tg s 



^0 



Here, ho/ 7k? is the integral 



h 
h G 



-x 2 /g s -x s /(3g s ) dx 



'2lX j _oo 

understood as a formal power series in g s . This can be explicitly calculated 



h G ttV2^ 3 2fc (2fc)! 



(B.19) 



(B.20) 



(B.21) 



(B.22) 



'B.23) 



(B.24) 



Notice that, as pointed out in [8] in the case of a quartic potential, the g s expansion of 



21 °g77J - logH(0,£ s 

rig 



(B.25) 



contains only odd powers of g s , as required in order to have an expansion of the free 
energy in terms of even powers of g s . Notice that the explicit solution for the F g in the 
one cut case provides an additional boundary condition for the holomorphic anomaly of 
matrix models . 
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C Modular forms and elliptic integrals 



We follow the conventions in [TT] . The complete elliptic integral of the first kind is defined 

as 

f 1 dt 

K(k)= / . (C.l) 

Jo v/(l-t 2 )(l-fc 2 * 2 ) 

The parameter k is called the elliptic modulus. Further one defines the complementary 

modulus as k' 2 = 1 — k 2 . The complete elliptic integral of the second kind is defined as 



"1 h _ 1,2+2 

dtJ T -^. (C.2) 



o 



The complete elliptic integrals of the first and second kind are related to each other by 
derivation, 

dK E{k)-k' 2 K{k) dE E{k) - K{k) 



dk kk' 2 ' dk k 

Useful transformation formulae are 

1 + *-*(*), 



(C.3) 




2 

^{E{k) + k'K{k)\ {CA) 
(l + k)K(k), 



as well as the Legendre relation 



E(k)K(k') + E(k')K(k) - K(k)K(k') = -. (C.5) 
Consider an elliptic geometry of the form 

4 

y 2 = l[( X - Xi ), (C.6) 

i=l 

where X\ < X2 < x% < are the branch cuts. Define the half-period ratio of the elliptic 
geometry, r, and the elliptic nome as q = e 17TT . It can be shown that 
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and moreover that 



n 02 7 2 _ ^2 r /2 _ ^4 



K(k) = -K e = ^ ^ = -f. (Ci 



2 d ' tf 4 



Here are the Jacobi theta-functions defined by 

$ 2 = Y^q^ n+1 * ] \ ^ 3 = ^g^ 2 , ^4 = ^(-l) n gK (C.9) 



The Thomae formula relates the branch cuts of an elliptic curve to theta-functions 
24"] . For the geometry consider above flC.6j) we obtain 



$2(7-) = -/C 2 0i - x 2 )(x 3 - x A ) 

$I(t) = -!C 2 (x 1 -x 4 )(x 2 -x 3 ) (CIO) 
•&\{r) = -/C 2 (xi - x 3 )(x 2 - x 4 ), 

and thus 

^A r )=^n^-^) 2 > (en) 

where K, and /C' are given in (I3.46p . 
It is convenient to introduce 

b = ^ 2 , c = #l d = #i, (C.12) 

where either two of them span the ring of T(2) modular forms. Here the congruence 
subgroup T(2) C SL(2, Z) is defined by 

r(2) = {7 G SL(2,Z) | 7 = 1 mod 2}. (C.13) 
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